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/ 'Using a simple model for outgassing from a small flat surface area, the sublimation rates of 

r carbon monoxide and carbon dioxide, two species more volatile than water ice that are known 
to be present in comets, are calculated for a suddenly activated discrete source on the rotating 
nucleus. The instantaneous sublimation rate depends upon the comet’s heliocentric distance 
and the Sun’s zenith angle at the location of the source. The values are derived for the constants 
of CO and C0 2 in an expression that yields the local rotation- averaged sublimation rate as 
a function of the comet’s spin parameters and the source’s cometocentric latitude. 



Numerous examples demonstrate that comets can flare up unpredictably at virtually any point 
of the orbit, including large heliocentric distances. An outburst of major proportions was recently 
experienced by Halley’s comet at 14 AU from the Sun (West et at. 1991). While the actual cause 
of these unexpected events cannot be stated with certainty, one of the more attractive amongst 
the plausible scenarios is sudden activation of a discrete source of dust ejecta (that makes up the 
observed halo) driven by expanding gases much more volatile than water (that have no transitions 
in the observed optical region of the spectrum and thus remain undetected). 

Judging from information based in part upon the results of experiments on board the space- 
craft that intercepted Halley’s comet, in part upon remote observations, carbon monoxide (e.g., 
Eberhardt et al. 1986, 1987, Woods et al. 1986, Krankowsky and Eberhardt 1990, Krankowsky 
1991) and carbon dioxide (Krankowsky et al. 1986, Feldman et al. 1986, Combes et al. 1988, 
Krankowsky and Eberhardt 1990, Krankowsky 1991) are the most likely species potentially capa- 
ble of instigating such eruptive events. 

To facilitate studies of the production of these volatiles from an isolated source on the nucleus 
surface, it is necessary to determine the diurnal and seasonal variations in the sublimation rate, 
that is, its dependence upon the Sun’s local zenith angle z and the comet’s heliocentric distance 
t. A simple model for outgassing from a flat surface element, based upon the assumption that the 
absorbed solar energy is spent on sublimation and thermal reradiation (but not on heat conduction 
into the nucleus) and already employed in studies of the water production (Sekanina 1988), has 
now been applied to carbon monoxide and carbon dioxide. The dependence of the saturated vapor 
pressure p (in dyn/cm 2 ) upon the temperature T (in °K) under equilibrium conditions has been 
expressed by 

log p = A-^ + CT, (1) 

where A = 10.11, B = 334, and C = 0 have been adopted for CO to fit the data listed by 
Egerton and Edmondson (1928), while A = 11.4320, B = 1275.62, C = 0.006833 for C0 2 have 
been taken from the same reference. The vapor pressures listed by Stull (1972) are generally in fair 
agreement with the adopted values of the coefficients, even though more elaborate formulas could 
be developed to improve the match. As in the case of water ice, solutions for the sublimation rate 
Z(z,r ) (measured in molecules/cm 2 /s) that was calculated from the equation of energy balance 
have been sought to satisfy a law 


Z(z,r)= Z 0 {r)-C{z,r), 


( 2 ) 
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where Z 0 is the sublimation rate at the subsolar point and £(z, r) < 1 is the dimensionless relative 
sublimation rate at the Sun’s zenith angle z , plotted for CO in Fig. 1 and for C0 2 in Fig. 2 at 
four different heliocentric distances. When almost all of the absorbed solar energy is spent on 
sublimation, the rate ( varies essentially as cos z. This is true up to a certain zenith distance of 
the Sun, beyond which ( begins to drop more rapidly. The approximation introduced in Sekanina 
(1988) for the relative sublimation rate of water ice has also proven useful for CO and C0 2 . Thus, 

£(z, r) = cos z — /(r) • sin 2 z for 0<z<z c , ( . 

= 0 for z>z c . ' ' 

Here z c = arccos{[l + (2/)" 2 ] 1 / 2 - (2/) -1 } is the Sun’s critical zenith angle, beyond which the 
sublimation rate is negligibly low compared with that at the subsolar point and /(r) is a function 
of heliocentric distance that is discussed below. The optimum representations of the relative 
sublimation rates that were achieved by applying the approximation (3) are shown in Figs. 1 and 2 
by dashed curves. This approach allows one to write the rotation- averaged sublimation rate ( Z ) 
at a given distance r (assumed constant during rotation) as 

(Z) = — / Z(z)dQ = ( Zo/tt ) -(a0 c + 6 sin0 c + c sin20 c ), (4) 

2t J~e c 

where 0 c (r) is a critical hour angle of the Sun that depends upon z c and the cometocentric latitudes 
of the active source, <£ ac t> and the subsolar point, 4 >ss? 

cos 0 C = cos z c sec 4>&ct sec <f> M — tan 4>&c t tan ^ 8S , (5) 



Fig. 1. Relative sublimation rate ((z,r) of carbon 
monoxide vs, the Sun’s zenith angle at the loca- 
tion of the source as a function of heliocentric dis- 
tance. The solid curves are the rates derived from 
the energy-balance equation. The broken curves are 
the least-squares solutions that accommodate the 
approximation (3). The curves have been shifted 
vertically for the sake of clarity. 



FIG. 2. Relative sublimation rate C( z > r ) of carbon 
dioxide vs. the Sun’s zenith angle at the location of 
the source as a function of heliocentric distance. 
The solid curves are the rates derived from the 
energy-balance equation. The broken curves are 
the least-squares solutions that accommodate the 
approximation (3). The curves have been shifted 
vertically for the sake of clarity. 
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FlG. 3. The sublimation rate 
per unit outgassing area at the 
subsolar point as a function of 
heliocentric distance for carbon 
monoxide and carbon dioxide, 
compared with the sublimation 
rate for water ice. 


and the coefficients a, 5, and c are functions of <^ act , </>ss? and f(r ): 

a = sin <£ act sin <p ss -\f{ 1 + sin 2 cos 2 <f> ss + cos 2<f> &ct sin 2 </> ss ), 

b = cos <?!> act cos <f> ss (l + 2/sin^actsin^gg), (6) 

c =\f COS 2 (/'act COS 2 (^ss- 

The peak (subsolar-point) sublimation rate Zo(r) is expressed by a slightly generalized version of 
the formula used in Sekanina (1988) for water ice: 

log Z 0 (r) = A + B log (r/r 0 ) + C (r/ro)^ + D log[l + (r/ r 0 ) 7 ] , (7) 

where A, 2?, C, D, /3, and 7 are constants for the given ice and tq is a “normalizing” heliocentric 
distance that is related to the heat of sublimation L (crudely, ro oc 1 /L 2 ). Figure 3 compares the 
heliocentric variations in the sublimation rates of CO and CO2 at the subsolar point with that of 
water ice. For CO the function /(r) is fitted empirically by 

f(r) = E (r/r 0 ) 2 [l + F(r/r 0 ) + G (r/r 0 ) 2 ] . (8) 

For CO2 this expression yields a fit that is somewhat inferior to that offered by a law 

/(r) = E (r/r 0 ) 7/4 exp [f (r/r 0 f] . (9) 

When t is in AU and Z$(r) in molecules/cm 2 /s, the coefficients have the numerical values listed 
in Table I for an assumed unit bolometric emisslvity and an albedo of 4 percent. The coefficients 
A through 7 provide an excellent fit up to r « 500 AU for CO and up to r « 40 AU for CO2. 
The coefficients E through G offer a very good fit up to r « 300 AU for CO and the coefficients 
E through G are applicable up to r « 20 AU for CO2. 


TABLE I. Sublimation constants for carbon monoxide and carbon dioxide. 

Species r 0 (AU) A BCDf3yEFGEFG 


CO 285 14.184 -2.00 -0.74 -0.63 2.09 14.95 1.386 -1.337 1.343 ... ... ... 

C0 2 20.2 15.961 -1.95 -0.75 -1.74 1.50 8.55 1.882 -1.775 1.779 0.78 0.90 2.70 
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Once the nucleus rotation parameters (determining the values of © c and cf> ss ) are available, 

Eqs. (2) through (9) describe fully the sublimation rate for any cometocentric latitude </> act , subject 
to the assumptions of the model employed. 

ACKNOWLEDGEMENT-' “ 

This research was carried out by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space Administration. 

;; REFERENCES : 

- ... _• I 

Combes M., Moroz V.I., Crovisier J., Encrenaz T., Bibring J.P., Grigoriev A.V., Sanko N.F., i 

Coron N., Crifo J. F., Gispert R., Bockelee-Morvan D., Nikolsky Yu. V., Krasnopolsky V. A., j 

Owen T., Emerich C., Lamarre J.M., and Rocard F. (1988) The 2.5-12 /xm spectrum of j 

comet Halley from the IKS-VEGA experiment. Icarus 76 , 404-436. { 

Eberhardt P., Krankowsky D., Schulte W., Dolder U., Lammerzahl P., Berthelier J. J., Woweries j 

J., Stubbemann U., Hodges R.R., Hoffman J.H., and flliano J.M. (1986) On the CO and 
N 2 abundance in comet Halley. In Exploration of Halley’s Comet, ESA SP-250 (B. Battrick, j 

E. J. Rolfe, and R. Reinhard, eds.), pp. 383-386. ESTEC, Noordwijk. j 

Eberhardt P., Krankowsky D., Schulte W., Dolder XL, Lammerzahl P., Berthelier J. J., Woweries I 

J., Stubbemann U., Hodges R. R., Hoffman J.H., and Hliano J.M. (1987) The CO and N 2 
abundance in comet P/Halley. Astron. Astrophys. 187 , 481-484. i 

Egerton A. C. and Edmondson W. (1928) Vapor pressure of chemical compounds in the crystalline 
state. In International Critical Tables (E. W. Washburn, ed.), vol. 3, pp. 207-209. McGraw- 
Hill, New York. s 

Feldman P.D., A’Hearn M.F., Festou M.C., McFadden L.A., Weaver H.A., and Woods T.N. 

(1986) Is CO 2 responsible for the outbursts of comet Halley? Nature 324, 433-436. 

Krankowsky D. (1991) The composition of comets. In Comets in the Post-Halley Era (R.L. 

j Newburn, Jr., M. Neugebauer, and J. Rahe, eds.), pp. 855-877. Kluwer, Dordrecht. 

Krankowsky D. and Eberhardt P. (1990) Evidence for the composition of ices in the nucleus of J 

comet Halley. In Comet Halleyi Investiyations , Results , Interpretations (J.W. Mason, ed.), 
vol. 1, pp. 273-289. Ellis Horwood, Chichester. 1 

Krankowsky D., Lammerzahl P., Herrwerth I., Woweries J., Eberhardt P., Dolder U., Herrmann 
U., Schulte W., Berthelier J. J., Uliano J. M., Hodges R. R., and Hoffman J. H. (1986) In situ 
gas and ion measurements at comet Halley. Nature 321 , 326-329. | 

Sekanina Z. (1988) Outgassing asymmetry of periodic comet Encke. I. Apparitions 1924-1984. ] 

Astron. J. 95 , 911-924, 970-971. 

Stull D.R. (1972) Vapor pressure. In American Institute of Physics Handbook (D.E. Gray, ed.), j 

p. 4-313. McGraw-Hill, New York. 

West R.M., Hainaut 0., and Smette A. (1991) Post-perihelion observations of P/Halley. III. An = 

outburst at r = 14.3 AU. Astron. Astrophys. 2 46, L77-L80. 

Woods TX, Feldman P.D., Dymond K.F., and Sahnow D.J. (1986) Rocket ultraviolet spec- 
troscopy of comet Halley and abundance of carbon monoxide and carbon. Nature 324, 436- 
438. 

1 1 

I ... :r ' -I... . - 1 

| | 


m 



